This document describes and presents initial results for the health-impact assessment related to adding more native plants to City of Denver. We aim to estimate the number of deaths that would be prevented by adding more green space, specifically native plants, to the City of Denver. We estimate the number of deaths prevented corresponding to three broad greening scenarios, which were informed by conversations with local stakeholders. The scenarios are:
For all scenarios, we measured green-space exposure by the normalized difference vegetation index (NDVI) measured by the Landsat-8 satellite at a spatial resolution of 30 meters squared. To determine the target NDVI value, we measured the NDVI of a 100% native zone in the Denver Botanic Gardens on several cloud-free spring and summer days. The NDVI value was about 0.5 in this area, so we assumed that, by adding more native plants to places in Denver, the NDVI of those areas would change to 0.5.
We first measured the mean NDVI of each census block group on July 4, 2021. The weighted block-group-level mean is presented here, weighted by the proportion each 30 square-meter pixel covers by the census block group. For example, if half of a pixel overlaps the block group, it receives a weight of 0.5 in the weighted average. We also removed bodies of water before measuring NDVI.
We excluded census block whose baseline NDVI was above the native-plant threshold:
We measured NDVI in a 200-foot buffer, a 100-foot buffer, and a 50-foot buffer around all bodies of water in Denver. We downloaded bodies of water from OpenStreetMap (code here. We defined residential exposure to these riparian areas as those individuals living within a 500-meter buffer, following the green-space literature that has defined green-space exposure based on residential proximity. We estimated the number of people in this buffer by multiplying the population density of the census block group by the intersecting area.
The below map depicts mean NDVI in the portions of census block groups that intersect a 200-foot buffer as well as those pieces that intersect the part of the 500 m buffer that would not be intervened upon, i.e., the part between the 200-foot buffer and the edge of the 500 m buffer.
We spoke with representatives at the local Office of Green Infrastructure, and they described three categories of initiatives–some planned, some aspirational–that could include native or adapted plants.
We obtained spatial data on parking lots from the City of Denver Open Data Catalog. About 9% of the city is covered with parking lots:
| area_mi2_prkng | den_area_mi2_total | percent_area_covered_by_parking |
|---|---|---|
| 13.56897 | 155 | 0.0875417 |
We similarly measured the NDVI on the parking lots as well as the NDVI in the census block groups that are within a 500-m buffer radius of any parking lots (which is most).
We present a small subset near Union Station, as the file size visualizing the NDVI measurements for the full city exceeds GitHub’s file-upload limits:
We estimated the number of deaths averted under each scenario by following a recent meta-analyses that estimated that for every 0.1 unit increase in exposure to NDVI, the relative risk of death decreases about 4% (pooled risk ratio of 0.96, 95% confidence interval (CI): 0.94, 0.97)1. We defined NDVI exposure as that within a 500-meter buffer of the individual’s place of residence. We used this risk ratio to estimate the population attributable fraction corresponding to the proposed change in NDVI from the baseline level to the native-plants level for each census block group (first scenario) or the part of the census block group intersecting the proposed intervention area (riparian and parking scenarios). We then estimated the number of deaths prevented considering estimates of population and all-cause mortality rates. We used 5-year American Community Survey Data (2015-2019) to estimate the population in 5-year age groups in each census block group and age-stratified estimates of mortality rates for the state of Colorado from the Global Burden of Disease project. If an intervention area did not cover an entire census block group, we multiplied the population density of the block group in that age group by the area of the block group covered by the intervention area to estimate the total population in the area affected by the intervention. We restricted analyses to adults aged 30 and above following the age range of many of the cohort studies reviewed. We also restricted analyses to those census block groups or census-block-group pieces below the native-plants NDVI threshold.
| scenario | scenario_sub | pop_affected | attrib_deaths |
|---|---|---|---|
| all-bg | 20-pct | 367077.00 | -44.0845890 |
| all-bg | 100-pct | 367077.00 | -227.6925261 |
| riparian | 200-ft | 171750.25 | -7.3772964 |
| riparian | 100-ft | 171750.25 | -2.4062585 |
| riparian | 50-ft | 171750.25 | -1.0475582 |
| ogi | ogi_proj-pct-75 | 72584.88 | -1.0589232 |
| ogi | parcel-pct-50 | 295353.72 | -0.7628711 |
| prkng | pct-prkng-100 | 360114.27 | -40.0231941 |
| prkng | pct-prkng-50 | 360114.27 | -19.9119495 |
| prkng | pct-prkng-20 | 360114.27 | -7.9410369 |
Add equity component
Add economic-impact attributable to health
Incorporate uncertainty throughout (via bootstrapping)
Copyright © 2022 Michael D. Garber